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I. INTRODUCTION 

The central object in the path-integral formulation of 
quantum statistics is the (Euclidean) transition ampli- 
tude A{a,b;t) = (b\e- t6 \a) 0, H S, 0. The starting 
point in setting up this formalism is the completeness 
relation 

A(a,b;t) = fdqi-- f dqN-i 

xA(a,q 1 ;s)---A(q N - 1 ,b;e), (1) 

where e = t/N denotes the time-slice width. To leading 
order in e the short-time transition amplitude reads in 
natural units 



A{q n ,q n+ x)£) 
1 



(27T£) d / : 



■ exp 



(gn+l - gn) 2 

2e 



eV(x n ) 



(2) 



where the potential V is evaluated at the mid-point co- 
ordinate x n — (q n + q n+ i)/2. Substituting @ in the 
completeness relation ([lj , the deviation of the obtained 
discrete amplitude from the continuum result turns out 
to be of the order O(e). This slow convergence to the 
continuum is the major reason for the low efficiency of 
the ubiquitous Path Integral Monte Carlo simulations [j| , 
especially in numeric studies of Bose-Einstein condensa- 
tion phenomena [H, 0, HI , quantum phase transitions and 
phase diagrams at low temperatures 0, Thus, in 

order to accelerate numerical calculations for statistical 
properties of quantum systems, it is indispensable to de- 
velop more efficient algorithms. The existence of such 



algorithms has been established recently [11 1 . 

To this end we worked out and numerically verified in 
a series of recent papers [H, [IH, Q, EEJ an efficient ana- 
lytical procedure for improving the convergence of path 
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integrals for single-particle transition amplitudes to the 
order 0(e p ) for arbitrary values of p. This was achieved 
by studying how discretizations of different coarseness 
are related to a hierarchy of effective discrete-time ac- 
tions which improve the convergence in a systematic way. 
In Ref. [H| we presented an equivalent approach which 
is based on a direct path-integral calculation of interme- 
diate time amplitudes to the order 0(e p ). The inherent 
simplicity of these direct calculations made it possible 
to extend the procedure to general many-body theories 
in arbitrary dimension and to obtain explicit results for 
the effective actions up to level p = 5. It turned out 
that increasing p leads to an exponential rise in com- 
plexity of the effective actions which, ultimately, limits 
the maximal value of p one can practically work with. 
These limitations of existing approaches, in particular 
in the case of many-body theories, are still below the 
calculational barrier stemming from this rise in complex- 
ity. This is a strong indication that new and more effi- 
cient calculational schemes must exist which should con- 
siderably improve the convergence of path integrals for 
general many-body theories. The availability of analytic 
expressions for higher p effective actions is essential for 
numerical calculations of path integrals with high pre- 
cision. Obtaining the information on energy spectra is 
just one important example of calculations that require 
high-precision numerical results. Furthermore, since the 
structure of higher order terms of effective actions is gov- 
erned by the quantum dynamics of the system, it can be 
used for extracting analytical information about the sys- 
tem properties. 

As is well known, in concrete calculations it is al- 
ways easier to solve the underlying Schrodinger equation 
than to directly evaluate the corresponding path inte- 
gral. For instance, in the case of particular potentials 
the Schrodinger equation approach allows an efficient re- 
cursive scheme to calculate perturbation series up to very 
high orders [13, [H, [H, H3, . With this in mind, we 
develop in the present paper a new and more efficient 
recursive approach for deriving the short-time transition 
amplitude from the underlying Schrodinger equation. To 
this end we proceed as follows. Section 2 presents the 
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instructive case of a single one-dimensional particle mov- 
ing in a general potential. The maximal level obtained 
by the new method is p = 35, and thus compares favor- 
ably with the best result p = 9 of previous approaches. 
In Section 3 we focus on the restricted problem of evalu- 
ating the velocity independent part of the discrete-time 
effective potential. We derive the differential equations 
for the velocity independent effective potential of a sin- 
gle one-dimensional particle and solve them for the case 
of a general potential up to level p — 37. Section 4 ex- 
tends the results of Section 2 to the case of a general 
many-body theory in d dimensions. We derive the dif- 
ferential equations for the general many-body effective 
potential and solve them up to level p = 10. Both the 
equations and their solution are presented in Section 5 
in a diagrammatic form which is similar to the recursive 
graphical construction of Feynman diagrams worked out 
in the series [H [H [H, [H, [H, [27[ . With this we illus- 
trate the inherent combinatoric nature of determining the 
discrete-time effective potential. The growing number of 
diagrams with level p leads to an increasing complexity 
of the expression for the discrete-time effective potential. 
We end by commenting on how this growth in complexity 
limits maximal attainable values of p. Throughout the 
paper we describe some important envisaged applications 
of the derived discretized effective actions for calculating 
statistical properties of relevant quantum systems. These 
applications include not only numerical calculations but 
also a new analytical approximation scheme which was 
made possible through the availability of high level effec- 
tive actions. 



where we have introduced V± = V{x ± x) as an abbrevi- 
ation. Their solution may be written in the form 



A = 
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: exp 
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where the effective potential W(x,x;e) is an even func- 
tion of x due to the symmetry ([3]) of the Euclidean transi- 
tion amplitude. Note that Eq. @ represents an approx- 
imation to the exact form §§§ up to order O(e). Substi- 
tuting © in and © yields 



W + x^+e^-'-ed^W-'-ed^W 
ox oe 8 8 

+ \ e 2 (dWf + I s 2 (dW) 2 = \ (V + + V.) , (10) 

xdW--eddW+-e 2 dWdW=-(V+-V). (11) 
4 4 2 



Both partial differential equations determine the effec- 
tive potential W(x,x;e) and thus the transition ampli- 
tude A(q, q';e). The initial condition © implies that W 
is regular in the vicinity of e = 0, i.e. it may be expanded 
in a power series in e. We are interested in using W to 
systematically speed up the convergence of discrete am- 
plitudes with N time slices to the continuum limit. This 
is done by evaluating W to higher powers in e. According 
to Eq. $2$ the dominant term for the short-time propa- 
gation is the diffusion relation x 2 cx e. Therefore, we 
expand W in a double power series in both e and x 2 : 



II. ONE PARTICLE IN ONE DIMENSION 

We start with calculating the transition amplitude 
A(q,q';e) for one particle in one dimension. It obeys 
the symmetry 

A(q,q'; S )=A(q',q;e) (3) 
and satisfies the time-dependent Schrodinger equations 



dS~2W 2+v{q) 



A{q,q';e) = 0. 
A(q,q';e) = 



with the initial condition 

A(q,q , ;0) = S(q-q'). 



(4) 
(5) 

(6) 



In terms of the deviation x = (q' — q)/2 and the mid-point 
coordinate x = (q + q')/2, both equations are rewritten 
according to 



A 



0, (7) 



[- dd + 2(V+-V-)] A 



0. 



W(x, x; e) 
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(12) 



Restricting the sum over m from to p — 1 leads to a dis- 
crete amplitude that converges to the continuum result 
as e p , i.e. as 1/N P . For later convenience we define all co- 
efficients c m .ki which are not explicitly used in Eq. \\2\ , 
to be zero, i.e. we set c m ,k — whenever the condition 
to > k > is not satisfied. 



Substituting the expansion (|12|) into the partial dif- 
ferential equations (|10[) and (fTTj) leads to two equivalent 
recursion relations. The second recursion relation turns 
out to be more difficult to solve to higher orders as it 
directly determines not the coefficients c m ,k{x) but their 
first derivatives. For this reason we restrict ourselves in 
the remainder of this section to the recursion relation fol- 
lowing from Eq. (fTO)) . The diagonal coefficients are given 
by 



y(2m) 

"m,m 77^ : 7TT j 
(2771 + 1)! 



(13) 



(8) while off-diagonal coefficients satisfy the recursion rela- 
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tion 



(m + k + 1) c m , fc = (2fc + 2)(2k + 1) c mjfc+1 

m-2 



~l~ C m— l,fc — C i-r c m-l-2,k-r 

1=0 r 

m-2 

- X! H 2 r ( 2fc ~ 2? " + 2 ) ^ Cm-i-l,fc-r+l > (14) 
(=1 r 

where the sum over r goes from max{0, k — m + 1 + 2} to 
min{fc, i} in accordance with the restriction that c m ^ = 
whenever the condition m > k > is not satisfied. For a 
given value of to, the coefficients c m ^ for k = 0, 1, . . . , to 
are determined as follows. The diagonal coefficient c mjI „ 
is directly given by (I13|) . whereas the off-diagonal coef- 
ficients c m ,k follow recursively from evaluating (|14|) for 
fc = to — 1, . . . , 1, 0. This recursive solution method is 
schematically depicted in Fig. [TJ Let us illustrate this 
procedure for the lowest levels. For p = 1 we immedi- 
ately obtain from (fl3"]) 



co.o = V . 



(15) 



For p — 2 we have to first determine ci,i from (|13l) . yield- 
ing 



Cll 



6 



(16) 



Then recursion relation (|14p states that c^o follows from 
the previously determined coefficients according to 
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From (fTS")) -(117 [) we then read off the result 
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Similarly, we find for p 
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The outlined procedure continues in the same way for 
higher levels p. We have automatized this procedure and 
implemented it using the Mathematica 6.0 package [28| 
for symbolic calculus. Using this we determined the ef- 
fective action for a one-dimensional particle in a general 
potential up to the level p = 35. Although the effective 
actions grow in complexity with level p, the Schrodinger 
equation method for calculating the discrete-time effec- 
tive actions turns out to be extremely efficient. The ul- 
timate value of p = 35 far surpasses the previously ob- 
tained best result of p = 9, and is only limited by the 
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FIG. 1: (Color online) Order in which the coefficients c m ,k are 
calculated. Diagonal ones follow from Eq. (|13[) , off-diagonal 
ones from the recursion relation (1141). 



sheer size of the expression for the effective action of a 
general theory at such a high level. 

The whole technique can be pushed much further when 
working on specific potential classes. For example, for a 
particle moving in a quartic potential we have obtained 
effective actions up to p = 140. Similar levels have been 
achieved for higher-order polynomial potentials. The in- 
crease in level and the decrease in size of the expressions 
for the effective actions originate from functional rela- 
tions between the potential and its derivatives. These 
relations are particularly simple in the case of polyno- 
mial interactions where all the derivatives of the poten- 
tial above a certain degree vanish. However, the benefits 
of working within a specific class of potentials are not 
only limited to polynomial interactions. For example, the 
functional relations for the modified Poschl- Teller poten- 
tial have allowed us to obtain effective actions up to level 

p = 40 dU. 

As already stated, the principal rationale behind con- 
structing high-level effective actions is to use them for 
speeding up Monte Carlo calculations. However, having 
obtained explicit expressions for effective actions to such 
high levels, it now becomes possible to use them exten- 
sively in both numerical and analytical calculations. In 
particular, having obtained an extremely precise knowl- 
edge of the short-time propagation of a system, it is pos- 
sible to use standard resummation techniques such as the 
Pade and the Borel method to extract useful information 
about its behavior for long propagation times. 

The derived effective actions can also be applied to 
systematically improve the Numerical Matrix Diagonal- 
ization (NMD) method [3(1 HH, [HJ for calculating energy 
eigenvalues and eigenstates. Note that the propagation 
time t used in the NMD method is just a parameter which 
is chosen in such a way that it minimizes the error associ- 
ated with the calculated energy eigenvalues. Therefore, it 
is always possible to select this parameter to be small, so 
that the obtained expansion of the ideal effective action 
can be used to substantially improve NMD calculations. 
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Furthermore, in this case we can even use analytic N = 1 
approximation for the path integral. In this approxima- 
tion there are no integrals to perform in Eq. ([1]) and the 
amplitude is directly given by the analytic expression ^ . 
Using such extremely rough discretizations is only pos- 
sible if one has determined the ideal effective action to 
very high orders p. In effect, one compensates without 
loss of precision the increase of discretization coarseness 
with the input of new analytical information concerning 
the propagation time which is contained in the effective 
action. In this way, without any integration or resumma- 
tion techniques, we can calculate large number of highly 
accurate energy eigenvalues, avoiding the usually needed 
limit t — ► oo which is difficult to approach. The large 
number of precise energy eigenvalues obtained by using 
this improved NMD method allows for calculating am- 
plitudes for longer propagation times using the spectral 
decomposition. This also makes it possible to calculate 
partition functions as well as global, local, and bilocal 
densities of states or other relevant statistical quantities 
with high accuracy. 



III. VELOCITY INDEPENDENT PART OF 
EFFECTIVE POTENTIAL 



The velocity independent part of the effective poten- 
tial Wo(x;e) = W(x,0;e) determines diagonal transition 
amplitudes. Although it does not contain all the informa- 
tion which is needed for constructing the effective action, 
it is of interest for determining physical quantities such 
as the particle density and the energy spectra. This ob- 
ject is also much simpler than the full W. In addition, 
the relation between V and Wo allows us to visualize the 
effects of quantization and discretization for a given po- 
tential. In this section we derive the differential equation 
for Wq for a single particle moving in one dimension. 

Both Eqs. ([7]) and ([8]) contain derivatives with respect 
to x, so it is impossible to set x — and obtain equations 
for Wo- However, differentiating Eq. (jHJ) with respect to 
x we get 



dd 2 A = AdVA 



(22) 



where the dots denote terms which vanish when x — ► 0. 
Thus, differentiating Eq. ([7]) with respect to x and sub- 
stituting the above result we obtain an equation which 
does not contain derivatives with respect to x. Finally, 
setting i = 0we find the differential equation for diago- 
nal transition amplitudes 



de 



1 



1 



dV 



Substituting A(x,x;e) 



Vd) A(x,x;e) 



(27re)- 1 / 2 exp(-eVF ) then 



. 



(23) 



yields the equation for Wq: 

dW + - ( 8 d^- - 8W dW Q + 8VdW - d 3 W ) (24) 
4 \ oe J 

-— 8dW —^- - 3dW d 2 W - — (dW ) 3 = dV . 
4 \ oe J 4 

This is solved in the form of the power series 

oo 

W (x;e) = J2cm,o(x)e m . (25) 
m=0 

Inserting this into the differential equation (|25[) deter- 
mines the coefficients c m ^o through the simple recursion 
relation 



(2m + 1) 4 >0 = i C-i, - 2 Vc' m _ h0 

■m — 1 771 — 1 



(26) 



+2 ^ c' k C m ^k~l,0 + 2 2^ ^ °k,0 c m-fe-l,0 
fc=0 k=l 
2 771—2 ^ m— 3 m—k—3 

~ J C fe,0 C m-fc-2,0 + T C k,Q C l,Q C r 



k=0 



^ /_ ^ £ ^ ! .' .U ""771 — k — I — 3,0 ' 

fe=0 1=0 



With this we have evaluated the velocity independent 
part of the effective potential up to level p = 37 for a 
particle moving in a generic potential V(x). As before, 
for specific potential classes one can go to much higher 
levels. 



IV. MANY-BODY SYSTEMS 

Now we extend the calculations of Section 2 to the 
case of a general non-relativistic theory of M particles 
in d dimensions. The derivation of the equations for W 
proceeds completely parallel to the case of one particle 
in one dimension. The Schrodinger equations now have 
the form 



M 



M 



A(q,q';e) = 0, (27) 
A(q,q';e) = 0, (28) 



where Aj and stand for c?-dimensional Laplacians over 
initial and final coordinates of particle i, while q and q' 
are dxM dimensional vectors representing positions of all 
particles at the initial and final moment. Furthermore, 
the potential V contains both the external potential and 
the respective interaction potentials between two and 
more particles. After substituting the rfM-dimensional 
generalization of the expression for the transition ampli- 
tude (JHJ) into the Schrodinger equations (|2"7) and (|28p , we 
find the analogues of Eqs. (fTU|) and (fTTj) for the effective 
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potential W 
W 



x ■ dW 
1 



dW 

de 



I ed 2 W - I ed 2 W 



+ ie 2 (dW) 2 + I e 2 {dW) 2 = \{V + + V-) , (29) 



x-dW - -d-dW- 
4 
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— (dW) ■ (dW) 



V) 



(30) 



Here we have used the definition A ■ B = AiBi, where 
2 = 1,2,..., Md and repeated indices are summed over. 
Either of the above two equations for W can be used 
to determine the appropriate short-time expansion as a 
double Taylor series in powers of e and even powers of x: 
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i d 2 W h0 + ^ d 2 W 2 ,i - i (dW ,of 



24 ' 24 

-i- d 4 V - i (dV) ■ (dV) , 
240 24 v ; v ; ' 

1 {x ■ dfv 
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W 3 ,i = ^d 2 W 2A + ±- B 2 W 3 , 2 



5040 
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T, d2 W %2 + ^B 2 W^ = ^{x.dfd 2 V„ 



40 



40 



-^(S-d)d i V(x-d)d i V, 



W{x,x;s) =^2^£ 

m=0 fc=0 



m — k 



W m , k (x,x) 



(31) 



x% x X{^ • ' ' Xi 2k c m (*£)• It turns 



where W m ^{x,x) 
out to be advantageous to use recursion relations for the 
fully contracted quantities W m< k rather than the respec- 
tive coefficients c^'' fc ' ,l2fe . In this way we avoid the com- 
putationally expensive symmetrization over all indices 
,?2fc- Again it is easier to work with the first of 
the two equations for W . Substituting (f3"Tj) into 
directly yields the diagonal coefficients 



W„ 



1 



(2m + 1)! 



{x ■ dy m V 



(32) 



The off-diagonal coefficients satisfy the recursion relation 
which represents a generalization of Eq. (fT4")l : 



8 (m + k + 1) W m , k = d 2 W m ^ k + d 2 W m , k+Y 

]{dW l<r ) ■ (dW m + 2 ,k-r) 



m-2 

ED 

1=0 

m-2 



J2 • w 



m—l—X,k—r-\-l 



)• 



(33) 



m 



As before, the sum over r goes from max{0, k 
I + 2} to min{fc,Z}. The above recursion disentangles, 
in complete analogy with the previously outlined case of 
one particle in one dimension. To illustrate this we write 
down and solve the equations up to level p = 4: 



W ,o 
Wi,i 

W 2 .2 

W 2il 



1 



G 



(x-dyv, 



120 ^ 9 ^< 

T2 d2w ^ + T2 B2w ^ = lk 



(x-d) 2 d 2 v, 



w 3fi = 



d 2 W 2fi + i a 2 W 3 ,i - (5Wq,o) • (SWi, ) 



360 J J 



V. GRAPHICAL REPRESENTATION 



(34) 



The above equations and their solutions can be cast 
in a diagrammatic form which is similar to the recursive 
graphical construction of Feynman diagrams worked out 
in the scries [H, H, H H [2(| . The effective poten- 
tial W represents the sum of all connected vacuum dia- 
grams of the underlying theory with the following Feyn- 
man rules. The propagator is represented by the Kro- 
necker delta 



'3 



(35) 



the Z-point vertex is the l-th derivative of the potential 



c\<9 i2 ■■■d il V = 



J*. 

12 



(36) 



and the even number of external sources stand for the 
discrete velocities x 



Xj, — 



(37) 



A simple dimensional analysis determines those dia- 
grams which contribute to a given coefficient W m .k- The 
discrete-time effective potential W is the generating func- 
tional of connected diagrams since it appears in the expo- 
nent. The diagrammatic notation makes it explicit that 
the short-time expansion of the effective discrete poten- 
tial is a purely combinatoric problem in which all the 
information is contained in the symmetry factors mul- 
tiplying individual diagrams. Thus, Eq. (I33|) represents 
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the Schwinger-Dyson equation of the underlying theory. 
As such it is the simplest way for actually determining 
the symmetry factors. 

The Schwinger-Dyson equation is now cast in a dia- 
grammatic form. To this end we begin with introducing 
general diagrams for W m ,k 



W„ 



m, k 



(38) 



2k 



According to (p?2"|) diagonal terms W m , m are directly given 
in terms of vertices which are contracted with even num- 
bers of external sources: 



W„ 



m, m 



^2m 



l 

(2m+l)! 



(39) 



2m 



The off-diagonal recursion relation (|33[) contains a differ- 
entiation of diagrams with respect to the mid-point co- 
ordinate x and the discrete velocity x. These operations 
act as follows: 



m, k 



m, k 
^2k 



,k 



2k 



(40) 



2 k 



m, k 
2k^- 1 



(41) 



Putting all those elements together we find the graphical 
representation of the Schwinger-Dyson equation: 



• + fc+l) 



m, k 



CL 

m — l,k 



+ (2fc + 2)(2fc + 1) 



0_ 

m, k + 1 



(42) 
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(=0 r 
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"~2r 



2k 

m — I — 2, k — r 



2 k 



2k^- 2r 



^^2r(2/c-2r + 2) 



Lr 



1=1 r 



2^- 1 



m — I — l,k — r + 1 



2k -2r+l 



The sum over r has the range as defined after Eq. (f3"3")) . 
From this we read off via complete induction that, in- 
deed, all vacuum diagrams contributing to the effective 
potential are connected. In this diagrammatic notation, 



W ,o 



Wl, = To 



W2.1 = rh 



f=|(i) 2 5 
# =lio(l) 2 (ll) I 

A 



j 



the previously obtained solutions for the discrete-time 
effective potential of a general many-body theory up to 
level p = 4 read: 



^2,0 - 230 



J_ 

24 



3,3 



W, 2 = 



3 > 2 3360 



240 (ll) 2 (12) , 

m =3M)(i) 4 (n) 
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5040 



1 

5040 
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IT: 



3,1 
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33(50 



07211 



VP 120 ^ flP 360 



^ (1) 2 (11) 2 -^(1)3(12)- i (1)(2)(12) 



33(50 



120 




1 

360 



i (11)3 i (H)(12) * (12) 



r 



Here we have also introduced a simple topological no- 
tation for the diagrams. The translation between both 
notations is obvious: (1) is an external source on vertex 
one, (11) a loop on that vertex, (12) is a link between 



vertices one and two. The topological notation makes it 
possible to present the effective action terms up to level 
p = 6 in a relatively compact form: 



— ^(l) 8 , 
362880 v ' 

— - — (1) 6 (H) , 
181440 ' ' 
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120960 ' ' 3360 w v ' 2520 ; 2016 w K ' v ' 

— - — (1) 2 (H) 3 — (1) 2 (11)(12) — (1) 2 (12) 2 — (1) 2 (22)(12) 

120960 v ' K ' 1680 v 2520 ^ ' v ' 1260 ^ 

— (1)(2)(11)(12) — (1)(2)(12) 2 , 

2520 1 A A A ' 5040 v A A ; ' 

W i0 = — - — (ll) 4 — (11) 2 (12) — (11)(12) 2 — (11)(22)(12) 

4,0 241920 V ; 2240 y ' K J 1680 1 A ; 40320 
1 i 



W 4 ,4 
^4,3 

W 4 ,i 



— — (12) 3 + — (12)(13), 
6720 v ; 240 v A J ' 



5:5 39916800 V ' 

W 5A = 

w 5 , 3 = 

Ws.2 = 



w 5A = 



(1) 8 (11) , 

15966720 w v ' 
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7983360 v ' v ' 181440 v ' K ' 90720 v 
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-90720 « W - 45560 WWDM - 45^60 W^ 12 )' 

-37i00 - Isk WW™ 12 ) - 50I0O 

5322^40 ^ W - iok W'* 11 ^) - 30240 C^ 11 * 12 ) 9 

(1) 2 (11)(22)(12) - ^ (1) 2 (22) 2 (12) - ^ (1) 2 (22)(12) 2 
1 „ 1 



~ 17280 ^ 

60180 WW - 60180 WX 11 )^ 2 ) " 302I0 M 2 )™ 12 )' 

-67200 W^ 11 ^ 12 ) - 10^00 W(2)(12)3 + sko 

~ (D 2 (12)(23) + JL (1)(2 ) (12 ) (13) + ^ (1)(2)(13)(23) , 
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' (12) 4 ~ 77^ (H)(22)(12) 2 - (H) 2 (22)(12) 



151200 ' ' 403200 ' ' v ' v ' 241920 
-2210 ("X 1 ^ 13 ) + ^60 f 11 )^)^) + i4 + 5^70 < 12 >( 23 >( 13 > 



Higher-level expressions are more cumbersome and may 
be found on our web site [2!|. Note that the diagram- 
matic notation reveals the fact that, as far as the short- 
time expansion is concerned, all systems fall into one of 
two classes of complexity depending on the value of the 
product dM. Discrete-time effective potentials for sys- 
tems with dM > 1 grow faster in complexity with in- 
creasing p than their dM = 1 analogues, as illustrated 
in Fig. [21 The reason for this is that several distinct 
diagrams collapse into a single term in the case of one 
particle in one dimension. Symbolical algebraic calcula- 
tions for dM > 1 effective actions were done using the 
program [291 ] written in Mathematica 6.0 in conjunction 
with the MathTensor package [33| . Using this program, 
for a general many-body theory we have derived analytic 
expressions for effective actions up to level p = 10. 

As in the case of one-dimensional systems, the obtained 
discretized effective actions can be applied to a host of 
relevant physical many-body problems in the framework 
of the Path Integral Monte Carlo a ppr oach, including the 
continuous-space worm algorithm [341 ] . For instance, our 
approach is applicable to efficiently determine the statis- 
tical properties of Bose-Einstein condensates which are 
confined in harmonic or anharmonic traps [HI, [36|, H3, HH ■ 



le+07 




1 5 10 15 20 25 30 35 40 

FIG. 2: (Color online) The number of diagrams contributing 
to a given level illustrates the rise in complexity with increas- 
ing p. Systems fall into one of two classes of complexity: the 
more complex ones with dM > 1 (top curve) and the less 
complex ones with dM = 1 (middle curve). In addition, the 
bottom curve gives the number of velocity independent dia- 
grams for dM = 1. 



Furthermore, our method is ideally suited for dealing 
with dilute quantum gases in a disorder environment 
where the impact of two-particle interactions upon the 
recently discovered phenomenon of Anderson localization 
is at present studied [39|, SfJ . 



VI. CONCLUSIONS 

We have given a detailed presentation of an analytic 
procedure for determining the short-time propagation of 
a general non-relativistic M-particle theory in d dimen- 
sions to extremely high orders. The procedure is based 
on recursively solving the Schrodinger equation for the 
transition amplitude in a power series of the propaga- 
tion time. This leads to a new recursion relation that 
has been solved to order p — 10 for the case of a general 
many-body theory. For a single particle moving in a gen- 
eral potential in d = 1 we have even achieved p = 35. In 
addition, for specific classes of potentials as, for instance, 
polynomial potentials, the equations have been solved to 
order p = 140. The resulting Schwinger-Dyson equation 
and its recursive solution have also been cast both in a 
familiar diagrammatic and a compact topological nota- 
tion. The presented results define the state-of-the-art for 
calculating short-time expansion amplitudes. They can 
be used to obtain orders of magnitude speedup in Path 
Integral Monte Carlo calculations. Thus, the extremely 
high orders of the short-time expansion make it possi- 
ble to perform new and precise analytical calculations 
of thermodynamical and dynamical properties of many- 
body systems. A list of applications of the presented 
method to relevant physical systems is briefly outlined. 
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